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Abstract 

A closed set of equations for the evolution of linear perturbations of homo- 
geneous, isotropic cosmological models can be obtained in various ways. The 
simplest approach is to assume a macroscopic equation of state, e.g. that of 
a perfect fluid. For a more refined description of the early universe, a mi- 
croscopic treatment is required. The purpose of this paper is to compare the 
approach based on classical kinetic theory to the more recent thermal-field- 
theory approach. It is shown that in the high-temperature limit the latter 
describes cosmological perturbations supported by collisionless, massless mat- 
ter, wherein it is equivalent to the kinetic theory approach. The dependence 
of the perturbations in a system of a collisionless gas and a perfect fluid on 
the initial data is discussed in some detail. All singular and regular solutions 
are found analytically. 
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I. INTRODUCTION 



Early progenitors of large-scale structure of the universe are usually provided for by small 
(linear) perturbations of otherwise homogeneous and isotropic cosmological models [l|] . After 
they have come within the Hubble horizon of the growing universe and their substratum has 
turned non-relativistic, they eventually become large, leading to gravitational collaps — the 
corresponding mass scale being set by the Jeans mass f2j. 

The theory of linear perturbations in a Friedmann-Lemaitre- Robertson- Walker (FLRW) 
model dates back to Lifshitz' work in 1946 ||. The metric perturbations and correspondingly 
the perturbed energy-momentum tensor are decomposed into scalar, vector and tensor parts 
according to their behavior under spatial coordinate transformations. In linear theory these 
parts evolve independently. Any tensor V* v is divided into a background part t^ v and a 
perturbation 5t^ v := t^ v — V L V . The background is given by a FLRW model. However, 
the mapping of points x^ on the physical manifold to points x^ on the background is not 
unique. Unphysical (gauge) modes can be avoided by using Bardeen's f|] gauge-invariant 
variables for the metric and matter components. The evolution of the metric and the matter 
perturbations is governed by the Einstein equation 

+ 5G^ = -8nG (f^ + 5T^ U ) , (1.1) 

where G^ v is the Einstein tensor and T fl u the energy-momentum tensor. Useful reviews of 
gauge-invariant cosmological perturbation theory are, e.g., Refs. and ||. 

In this paper we are mainly interested in comparing the methods of determining self- 
consistently the matter perturbations 5T fl u , needed to solve fll.ip. The Einstein equation 
describes the evolution of the metric perturbations for given matter perturbations. The 
covariant conservation of the perturbed energy-momentum tensor follows from the Bianchi 
identity. This does not fix 5T^ V completely; additional input is needed. The simplest 
possibility is to assume that both the perturbed and the unperturbed medium satisfy a 
simple macroscopic equation of state, e.g., that of a perfect fluid 

More complicated forms of matter in the early universe require of course a microscopic 
description. Such a description was first developed by Peebles and Yu || on the basis of 
kinetic theory || . There the Boltzmann equation in seven dimensional phase space (r, x, p) 
with an ansatz for the collision term is used, r denotes the conformal time, x the spatial 
coordinate and p the spatial momentum. The timelike component of the momentum is 
not independent, since the particles are on the mass-shell p° = \Jp 2 + m 2 (p = |p|). The 
Boltzmann equation together with Eq. ( |1 . 1|) has then usually to be solved numerically. A 
lot of work has been accomplished within this framework [|^|IT)|-|T5| . 

A different approach has been formulated more recently by Kraemmer and one of the 
present authors |16|]. In a completely field-theoretical framework, the matter perturbation 



bT^v is determined by the (thermal) graviton self-energy [fTTH , evaluated on the given back- 



ground. In contrast to kinetic theory, this framework is fully quantum theoretical from the 
beginning, and might turn out to be an alternative to the still to be formulated quantum 
kinetic theory. Up to now, this approach was used to incorporate collisionless, massless mat- 



ter Ijlgjlq] possibly mixed with a perfect fluid component p!9| , |20| , which can be described 
perturbatively by the high-temperature limit of the underlying field theory. In this limit 
one may expect equivalence with classical kinetic theory applied to the quanta of the field 
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theory. A detailed discussion of the relation between the thermal-field-theory approach and 
the one based on kinetic theory is however still missing and will be the main purpose of 
the present paper. The two approaches are compared for the case of collisionless, massless 
matter, with an arbitrary admixture of a perfect fluid. We will show that both approaches 
coincide in the high-temperature limit of thermal field theory. The exact solutions first found 
in the field-theoretical approach [[i6l , |T9| in the form of power series with infinite radius of 



convergence will be discussed in some detail with regard to their dependence on the initial 
data. Moreover, by a generalized power-series ansatz further exact solutions are obtained 
that are singular as the initial big-bang singularity is approached. (The existence of such 
solutions has been established first by Zakharov |I2[ and discussed further by Vishniac ||13|| .) 

The paper is organized as follows. In Sect. II we introduce the Einstein-de Sitter back- 
ground and the gauge-invariant metric potentials and matter perturbations. The Einstein 
equations are written down for these variables and their solutions for a perfect fluid universe 
are recapitulated. The kinetic theory is explained in Sect. III. We follow and add on to the 
approach of Kasai and Tomita [pl| , which defines gauge-invariant distribution functions. A 
set of integro-differential equations similar to that already existing in the literature |]10[ is 
derived for the most general initial conditions. In Sect. IV we prove that the same equations 
emerge from the thermal-field-theory approach in the high-temperature limit. The gauge- 
invariant equations describing a mixture of a collisionless gas and an isentropic perfect fluid 
are derived in Sect. V, and their exact solutions in terms of (generalized) power series (details 
are given in the Appendix) are discussed in Sect. VI. A summary is given in Sect. VII. 

Throughout the paper we use units % = c = ks = 1- Our notation is similar to that 
of H and |T5p, but we do not normalize energies and masses to v87rG. It differs only 



in the definition of the metric potentials for scalar perturbations and in that we expand all 
functions into planar waves instead of spherical harmonics, since we consider only a spatially 
flat FLRW model. Greek indices take their values in the set {0, 1,2,3} and Latin ones in 
{1,2,3}. 

II. GAUGE-INVARIANT COSMOLOGICAL PERTURBATIONS 

Before we define the gauge-invariant variables we introduce the Einstein-de Sitter back- 
ground. The invariant line element is given by 

ds 2 = S 2 {t) (-dr 2 + 5 ii dx i dx?) , (2.1) 

with r the conformal time measuring the size of the horizon in comoving coordinates Rh- 
The background energy-momentum tensor has perfect-fluid form, 

f» v = v?u v {E + P) + P6£ , (2.2) 

where u^u^ = — 1. We shall concentrate on the radiation dominated epoch, where E = 3P. 
The scale dependence of the mean energy density is given by E(S) = E(S = 1)S~ 4 , and the 
evolution of the cosmic scale factor S(t) follows from the Friedmann equation 



= ltoGE(S=l) T _ 
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As already explained in the introduction, the perturbed tensors are divided into scalar, 
vector and tensor parts: 



(2.4) 



Both tensors have ten independent components, four of which belong to the scalar, four to 
the vector and two to the tensor part (the latter correspond to the two polarizations of a 
gravitational wave). 

Since the Einstein-de Sitter background is spatially flat and depends only on the cosmic 
scale factor S(t), any function F(t, x) can be expanded into plane waves e lkx . A vector and 
a tensor valued plane wave is given by we ,kx and te lkx , respectively, where we specify the 
constant (but k-dependent) vector w to be transverse (w-k = 0) and the constant symmetric 
tensor t to be transverse-traceless (t • k = and tr t = 0). With this any symmetric tensor 
can be uniquely decomposed into its scalar, vector and tensor part. This is done for each 
individual mode proportional to e* kx of Eqs. ( |2.4| ) as follows [{/x, v} = {(0, i), (0, j)}]: 
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(2.6a) 
(2.6b) 
(2.6c) 



where k = k/fc, with k = |k|. 

The metric components A(r, k), B^ a \r, k), Hl(t, k), (r, k) and the matter variables 
8(t, k), v^(r,k), ttl(t, k) : 7iy (r,k) with a = {0,1} transform non-trivially under gauge 
transformations, i.e., coordinate transformations of the perturbed manifold while keeping the 
background fixed. Only the tensor components (a = 2) are gauge-invariant by themselves. 
The gauge transformations for scalar modes are 



r = r + T(r, k)e 



;kx 



x i — L^(t, k)ik l 



(2.7a) 



and 
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x i + L (1) (r, A;)wV kx 



(2.7b) 



for the vector perturbations. From ( 2.7a| ) and (|2.7b|) follows that only two scalar "potentials" 
and one vector "potential" of the metric perturbations are gauge-invariant. A set of such 
metric potentials is given by 



n 



h l + \hP + 1{b^-h t 



- A + H L 
for the scalar part and by 



(0) 



/(0) 



'(Oh 



T 



1{B^-H[ 



* := B« - Hp 



for the vector perturbations. The prime denotes a derivative with respect to 

rk. 



x 



(2.8a) 



(2.8b) 



(2.9) 



The variable x is related to the number of (half-) wavelengths inside the Hubble horizon 
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x 
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A 



R 



H, 



and we shall later adopt it as a normalized time variable. 
$ and II are translated into the variables used by Bardeen 



by 



<!> H 



2 



-n- 

2 



In the same manner gauge-invariant matter variables can be defined. For the scalar pertur- 
bations we introduce: 



7/ 



n L -5 
v<® - H. 



5 + ! 



/(0) 
T 

7j(°) - 



fi (0) 



- H 0) 



(2.10a) 



?7 is called the entropy perturbation, although it may not coincide with the perturbation in 
the true physical entropy, is the amplitude of the matter velocity, and is related to the 
shear a by a = v^k/S [§J. e m is the density contrast on a spacelike hypersurface which 
represents the local rest-frame of matter everywhere, whereelse e g is the density perturbation 
on a hypersurface whose normal vector has no shear. The anisotropic pressure tt^ is gauge- 
invariant by itself. 

For vector perturbations the invariant matter variables are 



,(i) - 



H. 



'(i) 



=: v c + * 



;2.10b) 



where Vf' can be understood as in the scalar case and v c is the velocity relative to the 
normal to the constant-time hypersurface. v c k/S is the intrinsic angular velocity (vorticity) 
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|| . The anisotropic pressure ir^ is gauge- invariant as above. For tensor perturbations the 
only physical variables are given by the metric potential Hj, =: H and the anisotropic 
pressure both already gauge-invariant. (Here we follow the conventions of Ref. [|J, 
which differ from the definitions of \I/ and H in Ref. |Tj| by a factor of 2.) 

From the Einstein equations (|1 . 1|) relations between the metric potentials $, II, and H 
and the matter variables T],e m ,e g ,ir^ ,v^ a \ and v c follow. The details are given in Ref. [|J. 
For the scalar components 

x 2 

y$ = e m (2.11a) 

x 2 n = 4 0) , (2.Hb) 

which shows the reason for our choice of the metric potentials: $ is a potential for the 
density perturbation and II the potential for the anisotropic pressure perturbation. Another 
set of equations is useful: 

x 2 ( $" + + + -IT - -u) = -77 (2.12a) 
V x 3 x 3 / 



(x 2 + 3)$ + 3x& + 6n = 3e g . (2.12b) 
Covariant conservation of the perturbed energy-momentum tensor implies 

^ (0) + = \ (e m + V ~ ~ \* ~ n- (2.13) 
For vector perturbations the Einstein equations read 

q = v c (2.14) 



,r 2 
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and conservation of the energy-momentum tensor yields 



< = -l4 ] - (2.15) 



Tensor perturbations fulfill 

X 2 (H" + -H' + H) =4 2) . (2.16) 



x 

Our task in the next Sec. (kinetic theory) and in Sec. IV (thermal field theory) is to 
self-consistently determine the matter perturbations from a microscopic theory. But the 
simplest case is to close the above equations by assuming macroscopic equations of state. 
We end this Section by briefly recapitulating the solutions obtained with a perfect fluid. 

In the case of a perfect fluid the anisotropic pressure vanishes (71^ = for all a). 



Additionally we assume adiabatic (isentropic) perturbations, i.e., rj = 0. Eq. (2.1 1] j implies 



II = 0, and Eq. (|2.12a ) reduces to an ordinary differential equation for the metric potential 
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$. The solution is given by a spherical Bessel function for regular behavior at x = and a 
spherical von Neumann function for singular behavior in the origin: 



where and are constants. The regular solution approaches a constant for small 
x (superhorizon scales), which leads to a growing e m ~ x 2 because of Eq. (|2.11| ). Vector 
perturbations in a perfect fluid have constant v c from ( [2.15| ), which is a consequence of the 
Helmholz-Kelvin circulation theorem |[21|| . Regular behavior of ^ at x = however implies 
w c = 0, for otherwise, \1/ would be proportional to 1/x 2 . The tensor perturbations are freely 
propagating gravity waves; the solution of ( [2.16| ) is given by: 



H 



c {2) Jo(x) 



yo{x) 



which can be regular or singular as x —>■ 0. Their amplitude decays 
accordance with energy conservation [plj . 



~ x 



for large x in 



III. KINETIC THEORY APPROACH 

The framework for kinetic theory within general relativity was set up long ago, see, 
e.g., Refs. jtf^l; a recent book on this subject is Ref. |23fl . In the context of cosmological 
perturbations kinetic theory was first used by Peebels and Yu ||. They discussed the coupled 
Boltzmann and Einstein equations for photons interacting with electrons via Thomson- 
scattering. For collisionless, massless neutrinos Stewart Jl0[ derived the full set of Einstein- 
Vlasov equations. An extensive investigation of density perturbations with massive and 
massless neutrinos, radiation and other matter was performed numerically by Bond and 
Szalay 



Kinetic theory is applicable for particles whose de Broglie wavelength is smaller than 
their mean free path. For collisionless matter this is certainly the case. Moreover, the de 
Broglie wavelength 2n/p has to be smaller than the Hubble horizon Rh = t ~ (GT 4 ) -1 / 2 , 
i.e., with 

p- 1 ~ T- 1 < m PL T- 2 , (3.1) 



the temperature T has to be below the Planck scale 

In this Section we shall concentrate on the case of collisionless, massless matter; the 
results for the more general and more realistic case of a two-component system including a 
perfect fluid will be given later on. The massless case is somewhat exceptional in that it 
is rigorously consistent with a thermal equilibrium situation. A thermal equilibrium in an 
expanding universe can only be obtained in the case of massless particles, since the FLRW 
models provide a conformal timelike Killing- vector only |§, or in the nonrelativistic limit. 

In setting up the kinetic approach to cosmological perturbations, we follow the paper 
by Kasai and Tomita Ul5 |. They use gauge-invariant distribution functions ||24j| , split up 



into scalar, vector and tensor modes. Other authors are usually working in a specific gauge 
T0|-14|. An infinite set of gauge-invariant equations for the fluid variables describing density 
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perturbations has been obtained in Ref. [p5|j . The authors of Ref. [D| use kinetic theory 



within the tetrad description of general relativity [g^j . The tetrad is defined by 
It spans an orthonormal basis in the tangent space. 

The phase space of an ensemble of particles is described by the four space-time coordi- 
nates x M and the momentum components of the particles p M := dx^/ds, s being an affine 
parameter. Since the particles' energy is on the mass-shell p° = |p| =: p (for m — 0), the 
dimension of the phase space is seven. An invariant volume of momentum 3-space is given by 
dp 1 dp 2 dp 3 /p =: d 3 p/p, where the indices refer to the tetrad basis. The distribution function 
F(t, x, p) fulfills the Boltzmann equation 

HF) ■= ~ VP P ^)j F ^ x, p) = C[F) . (3.2) 

The collision term C is a functional of F and has usually to be put in by hand. It is 
identically zero for collisionless matter. T^ p denotes the usual Christoffel symbol. 
The energy-momentum tensor is defined by 



T\ := / ^pPpuF . (3.3) 



In the collisionless case covariant conservation of T^ v follows from the Vlasov equation 
L(F) = 0. The perturbations of the energy-momentum tensor are given by 



8T» V = J^j- (p*PJF + {5p»p v + p^p v )F) . (3.4) 



The perturbations of the particle momenta 5p^ are calculated with help of the inverse tetrad, 
which reads 

1 1 



in linear theory. Sp 11 = 5e£p a relates 5p^ to the metric perturbations. In (|3.4|) the perturbed 
distribution function is 5F = F — F, where F = F(Sp) is the background distribution 
function. In thermal equilibrium, the latter is a Bose-Einstein or a Fermi-Dirac distribution 
depending on S(r)p, the energy of the particles in comoving coordinates, only. Thus 

1 



e ps P ± 1 



with (3S(t) being the inverse temperature. 

The decomposition into scalar, vector and tensor perturbations in the plane wave expan- 
sion is given by 



5F(r,x, P )=Ef/ (0) + / (1) ^ + / (2) 




(3.5) 



With an isotropic background distribution, the depend on r, the absolute values k 
and p, and the cosine \x := k • p/kp. For a given mode with wave vector k, we write = 
f( a \x,p, //). Like the matter variables introduced in Eq. ( p.6|) the distribution function is 
not invariant under the gauge transformations (|2.7f) . 
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A. Scalar perturbations 



For scalar perturbations a gauge-invariant distribution function is given by 



== / (0) +P% - Hp) + % h{\b® - Hp)) , (3.6) 



or by 



J{x,p,ri:=I®-^*. (3.7) 
According to the expansion theorem for Legendre functions, 

oo 

I<°W.aO = Y,a { °\x,p)i n PM (3.8) 

n=0 



for every E L 2 [— 1, 1]. The utility of ( |3.8| ) becomes apparent by calculating the gauge 

invariant matter variables e g , and tt^ from ( |2.10a| ) and ( |3.4j ) (a more detailed derivation 
is given in |15| ): 



4?r r°° - .3 jo). 



e a = -^ dpp ay(x,p) (3.9a) 
E Jo 



M = _l rdpp 3 a?(x,p) (3.9b) 
E Jo 



*t = T§ rdpp^\x,p) . (3.9c) 
5E Jo 

Here the orthogonality of the Legendre polynomials has been used. e m follows from e g and 
v(°\ 7] vanishes, since p^p^ = in Q3.4Q and S(p fl p fl ) = 0, as well. Inserting J into ( p.2|) for 
scalar perturbations leads to the equation 

1 dJ OF 
f ~ -P a ir-a + w J = -w + n ) , (3-10) 

where a takes the values 1, 2, 3 and refers to the tetrad basis. Its solution is 

r) TP rx 

J(x,p,(i) = e-^ {x - Xo) J(x ,p,fi) - tfip— / dx'($ + IT) (x')e-^ {x - x,) . (3.11) 

Op Jx 



With 



e g = i /— p 2 J{x,p,fi) -2$ 
E J p 



and 
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the Einstein equations (|2.12|) read: 



$" + -$' + $ = _($ + n)--($ + n)' (3.12a) 

x 3 x 



{x 2 + 3)$ + 3x$' = 6 ($ + II) - 12 [ X dx'j (x - x') ($ + n)' (x) + 12 f] /?n } j 

J 1(1 n 

(3.12b) 



x — xo J • 

X 



The spherical Bessel functions j n (x — xo) are the Fourier transforms of the Legendre poly- 
nomials in Eq. (|3.8|) . The last term on the r.h.s. determines the initial conditions. The 
coefficients j3^' are defined by 



4 0) := i/ >A?(%p)-nW-|$W 

/3(°) := f/ oo dp^oip)(x 0l p) , n>l. 



(i) 24tt 
5£ 



For vector perturbations Eq. (|3.2|) reads: 



(3.13) 



B. Vector perturbations 

Here we use the gauge-invariant distribution function 

and expand it, as above, into Legendre polynomials in \x 

oo 

J«(x,p,/i) = '£a<p(x, P )t n PM . (3.15) 

n=0 

In terms of the coefficients a^ 1 ) , the matter variables are 

v c = j^J o dpp 3 (a { o\x,p) + ^4 1} (x,p)^ , (3.16a) 



| / o °°rfpp 3 (-afW) + ^(^P)) • (3.16b) 



I'd) _ + = W— * • (3-17) 



Its solution is: 
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J (1) (x,p,/i) = e- l ^ x - x ^I^(x ,p,a)+iap— / dx'^ix')^ 1 ^-^ . (3.18) 

op 



■I'D 



With this solution inserted into v c , the Einstein eq. ( |2.14j ) becomes 

x 2 ^ = -2A [ X dx' j2{x - P- V(x') + 12 f; /#> ( Jn (x - x ) + %{x - s )) • (3.19) 

Jx XX , ;t=0 



The initial conditions for the Volterra type integral equation ( |3.19|) are given by the coeffi- 
cients 

& ] = - dpp 3 a£\x ,p) . (3.20) 
E Jo 



C. Tensor perturbations 

For tensor perturbations the distribution function is gauge-invariant itself. With its 
expansion into Legendre polynomials 

oo 

f^ 2 \x,p^) = J2 a£\x,pyPM (3.21) 

n=0 



the anisotropic pressure reads: 

i J^dpp 3 (a ( 2 \x,p) + 2 -a (2 \x,p) + I-af (x,p)) . (3.22) 



( 2 ) 24vr 

71 j. 



15E 

The Vlasov equation for the tensor perturbations reads 



x op a op 



(3.23) 



It is solved by 



f) TP rx 

f i2) (x, p, n) = e -*rt*-*o)_f (2) f \ + / dx'H\x')e-^ {x - x,) . (3.24) 

Inserting this into ( |2.16| ) yields 

x i ( H " + —H' + h) = -24 f dx' 3 -^—Q- H'(x') 

\ X J Jx {X — X'Y 

OO 

+3 E /$° - *o) + 2#(x - x ) + £ y) (x - x )) . (3.25) 

n=0 

The coefficients encoding the initial conditions are 



P$* = i r~dpp 3 aW(x ,p) . (3.26) 
E Jo 



oo 
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IV. THERMAL-FIELD-THEORY APPROACH 



Apart from some notable exceptions [^,^,^] all results given in the literature within 
kinetic theory have been derived by choosing a specific gauge. As was shown in the preceding 
section, this dependence on the gauge may be circumvented by completing the program lined 



out by Kasai and Tomita ||15|| . As we have seen, this requires a skillful redefinition of the 
basic distribution function. 

A manifestly gauge-invariant approach is provided by thermal field theory. For a system 
containing gravity as well as matter at a temperature T below the Planck scale (T <C mpi), 
the effective action (i.e. including all radiative contributions) T[g] can be split in a part 
describing classical gravity, the Einstein-Hilbert action S G itself, and an effective action 
r[g] A/ induced by the thermal matter with classical action S M , 

r\ T « mpL = s G + r M . (4.i) 

This is automatically gauge-invariant provided any gauge fields are subject to background- 
covariant gauge conditions. 

T A/ does not depend on curvature effects in the leading-order temperature contribution. 
In the radiation dominated regime, the Ricci tensor is proportional to GT A . Thus for temper- 
atures below the Planck scale (T <C itlpl) curvature effects are lower order in temperature, 
i.e. <T 2 . 

The energy-momentum tensor is given by 

2 ST M 

T^(x) := -;=-=— , (4.2) 
and self-consistent perturbations thereof have to fulfill 

(4 ' 3) 

This relates the perturbed matter variables to the thermal graviton self-energy which is 
defined by 

A perturbative expansion in Feynman diagrams is appropriate for weakly interacting 
matter. Higher loop orders due to internal graviton propagators are suppressed by powers 
of GT 2 <C 1 for T rripL. A high-temperature expansion of T M in the sense of k <C T is 
appropriate for the study of cosmological perturbations, because the external scale is set by 
the Hubble horizon Rh = t ~ (GT 4 ) -1 / 2 , which is 3> T _1 for temperatures well below the 
Planck scale. 

First attempts to calculate the thermal gravity self-energy have been undertaken in 
Refs. A complete calculation of the leading high-temperature contribution was given 
first by one of the present authors in Ref. JT7 ]. For this only one- loop diagrams without 



internal graviton lines need to be computed. Similar issues have been studied in Refs. [28 
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Recently, de Almeida et al. |29| have computed the next-to-leading order contributions for 



radiation and bosonic matter in the massive, but still collisionless case. 

By calculating the graviton self-energy from thermal field theory on the cosmological 
background, one can obtain the r.h.s. of the Einstein equations ( |1 . 1|) from ( |4.3|) using the 



definition ( |4.4j ). In the kinetic theory approach the Boltzmann-Einstein system of equations 
has to be solved simultaneously. This can be done analytically for collisionless, massless 
matter, but in general only numerical solutions can be obtained. In the thermal field theory 
approach, the corresponding problem is to calculate the graviton self-energy on a curved 
background space-time. An explicit reference to perturbations of distribution functions is 
completely obviated, which might turn out to be useful when more than the classical limit 
is of interest. A complete evaluation of the thermal graviton self-energy in curved space 
is certainly a formidable task. However, restricting our attention to collisionless, massless 
matter again, one can take advantage of the fact that the high-temperature limit of the 
effective action possesses, in addition to diffeomorphism invariance, an invariance under 
conformal transformations, which is expressed by the Ward identity |T6|JT9| 
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Since the Einstein-de Sitter background is conformally flat, the graviton self-energy can be 
simply calculated by evaluating 11^,5 on flat space-time and transferring it to the curved 
space by multiplication with cosmic scale factors S. This calculation was done by one of the 
authors in Refs. |i7|JT9|] . We will not repeat it here, but sketch the most important steps of 



its derivation. 

The graviton self-energy IIq,^ is related to the Fourier transformed self-energy in flat 
space-time Tlap-ys by 



U^(x,x')\ g=s2v = S~ 2 (t) J ^-/ Kx ~ x,) W vp \k) S~\r') . (4.5) 



v 



In the high-temperature limit, its tensorial structure turns out to be the same for any field 
theory and is given by 



= \^ vpu - ~ (6$Tf» + rJgTf*) Vysl^ 5 , (4.6) 



which follows in a straightforward calculation from the Feynman rules of the imaginary time 
formulation (see e.g. |3(|) of thermal field theory. The totally symmetric quantity l a ^ & is 
given by 

^.rEj^rtW. ( , 7) 

Sum and integral are conveniently performed with help of the so-called Dzyaloshinski algo- 



rithm |31|j30| . Its values are tabulated in the appendix of Ref. ^L7j. They are polynomials 
in k /k times 

n Ax 1 , fk + k\ 

g ° ( T ) = 2 ln U^ ' 
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the zeroth Legendre function of second kind. These may be rewritten with help of the 
formula 



n j 

Q n ( Z ) = P n (z)Q (z) - J2 -Pm-l(z)Pn-m(z) 



m=i m 



in terms of higher Legendre functions of second kind plus polynomials. 

Like all other quantities, the graviton self-energy is expanded into plane waves. Therefore, 
only the Fourier transformation of ko into r in ( |4.5|) has to be performed. The integrals 
involved are 

lim — / due-^ (x - x )Q n (cu) = (-i) n+1 e(x - x')j n (x - x') , 

7— >0+ ZlT J -00+17 

where uj = k$/k and retarded boundary conditions are imposed. 
In Eq. ( f4.3|) this leads to convolution integrals of the form 

/oo 
dx'Q{x - x')j n (x - x'){<$>, IT, H}'(x'). 

There has to be a finite non-negative integration bound xq because r > 0. On the other 
hand, the functional derivative in equation (|4.3|) is usually defined without restriction. For 
that reason we have to put in the initial conditions for xq by hand. We do this by the 
replacement of the metric potentials by 

00 

{$, n, vd, h}\x) - {$, n, H}'{x)e(x - Xo ) + E i^^H* - x Q ) . 

n=0 

The infinite sum provides the initial conditions for an formally infinite-order differential 
equation, which is obtained differentiating the convolution integral above (see |[19[). The 
infinite sum yields derivatives of spherical Bessel functions, which can be rewritten into 
spherical Bessel functions. 



A. Scalar perturbations 

To arrive at the expressions for the matter perturbations in thermal field theory, we have 
to relate them to the graviton self-energy through Eq. ( |4.3| ). The matter variables on the 
r.h.s. of Eqs. (|2.12|) read: 



r] = 



and 



2$ + 411 - 4 r dx'j (x - x) ($ + n)' (x') -4]T 7i 0) jJ n) ' 

Jx ° n=0 



x - X J 



(4.8) 



Inserting this into ( |2.12| ), the same equations as in the kinetic approach, Eq. ( p,12[ ), are 
obtained, except that the initial conditions are parametrized differently. The coefficients 
7^°' ) can be related to the (3^ appearing in Eq. (|3.12|) through the formula 
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' (ry'n-i - (n + l)j n+1 ) , (4.9) 



2n + 1 



by a linear transformation which is non-singular and of upper-triangular form (albeit infinite- 
dimensional) . If (3^ = for all n > 3, which is a frequently adopted simplifying assumption 



in the kinetic approach (the so-called 14-moment approximation, see, e.g., Ref. |3~2fl), then 
this relationship is given by 

7 f = -^f . (4.10) 



B. Vector perturbations 

In a similar manner we obtain for vector perturbations 



-3 / dx 



,J2{X-X 



7^(x') + 7i 1} Uo(x - x ) + j 2 (x - x )) 



(n) 



(4.11) 



n=0 



which together with ( |2.14| ) reproduces (|3.19|) , if the coefficients are re-expressed in 
terms of the (5^ which parametrize the initial conditions in Eq. ( 3.19Q . Because of jo + Jo = 
| (jo + J2), the relationship between the and the (5^> is, apart from an over-all sign, the 
same as in the scalar case. Hence, we have 



7$" 
7?' 



4" + \®\ 



(4.12) 



C. Tensor perturbations 



In this case 



7T 



(2) 



-24 



x'o (x — x') z 



n=0 



15" 



2 J 2 + 7^j4) (n) (^-^o) 
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35" 



(4.13) 



which with ( [2.1 6|) leads again toEq. (|3~25| ). Because of j + 2j '+j 7y) 
the first three coefficients and (3^ are related by 



7f 

7?' = « 
7f 



(2) 

3 
2 



,^Jo + ^j2 + ^j4), 



(4.14) 
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if all higher coefficients are zero, exactly as in the scalar case. 

Summarizing, the kinetic theory approach in the gauge invariant formulation of Ref. [15 
and the thermal-field-theory approach in the high temperatur limit lead to the same equa- 
tions for cosmological perturbations of a Einstein-de Sitter universe filled with collisionless, 
massless particles. There is a difference, however, in the way the initial conditions are intro- 
duced, which may seem somewhat arbitrary in the field theoretical approach. In fact, the 
initial condition terms in Eqs. (|3.12p , ( p.l9| ), and ( |3.25| ) are completely arbitrary functions 
due to the von Neumann expansion theorem. We could therefore just have added the in- 
homogeneous terms by hand at the end of the calculation, but as we have seen the way in 
which we have implemented them there is still a close relation to the inital moments of the 
initial distribution function of the kinetic approach. 



V. COLLISIONLESS MATTER AND A PERFECT FLUID 

The universe contains, besides collisionless particles, various different forms of matter as 
radiation and baryons. An example for a collisionless, massless gas is provided by neutri- 
nos after their decoupling below the electroweak scale. Another example are background 
gravitons, left over from the epoch of quantum gravity. In both examples all other matter 
may be described by a perfect fluid during a certain epoch, while collisions keep it close to 
thermal equilibrium. 

To model such systems we derive the equations for cosmological perturbations evolving 
in a mixture of a collisionless, massless gas and a perfect fluid. We assume that the gas and 
the fluid interact only through gravitational forces, thus 

T» v = T^ v + V„ . (5.1) 

To describe the radiation-dominated epoch of the universe, we assume that the components 
of the background perfect-fluid energy-momentum tensor T v ^ u obey the same equation of 
state as the collisionless gas components, i.e., E p{ = 3P p f. For both perturbed energy- 
momentum tensors, 8T C \\^ V and 8T P ^ V1 gauge- invariant matter variables are defined as in 
Sec. II. We define the ratio of the collisionless background energy- density to the total back- 
ground energy- density to be 

a := _ Ec \ . (5.2) 
E c \\ + E-pt 

In order to gain a closed set of gauge-invariant equations we allow only for adiabatic (isen- 
tropic) perturbations, i.e., r\ = ar) c n + (1 — a)r] pi = 0. Entropy perturbations would arise, 
e.g., from phase transitions ||, which lie beyond the scope of this work. 

The extension of Bardeens H] gauge-invariant formalism to a multi-fluid universe was 



given in Refs. [§,0]. In connection with a collisionless gas this was done in Ref. | 15fl . Since 
we deal only with two components, fulfilling the same background equation of state, the more 
general formalism of these references is not necessary here. The equations we are looking 



for have been obtained by one of the present authors in Ref. [ I9| , but without including the 
most general initial conditions. The following gauge-invariant equations are derived from 

(OD with (0). 
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A. Scalar perturbations 



To arrive at a closed set of equations for scalar perturbations we define $ p f + $ c u := $ 
with help of ( |2.11a| ), i.e., 



x 



^cii — o.e m en . 



(5.3) 



Eq. fl2.12a| ) does not change since we deal with adiabatic perturbations only. With ( ^.12a| ), 
(p.3|), and ( |2.11b| ), $, $ c n, and II can be calculated from e m c n and 7r^ dl . The anisotropic 
pressure of the perfect fluid is zero from its definition. Within the field-theoretic approach 
from Eqs. fl5J) and ( gTTJ ) 



x 



ell 



-2a 



$ - 2n + 2 r da;' (j (x - x') + -^(z - x')) ($ + n)'(x') 

Jxo \ X / 



+2^ 7 ( )( j0 + -j 1 ) (n) ^-n 1 ; 



n=0 



.r 



(5.4a) 



and 



x 2 n = -12a 



,(°)^)( x -a;o) 



dx'j 2 (a;-x')(<I» + n)'(a;') + E7i 0) jt 

-0 n=0 

follow. In terms of /3^'s (cf. Sec. Ill) the inhomogeneous terms respectively read: 

oo q 
4«£^ 0) (j«--Jn)(*-*0) 
n=0 X 

and 

oo 

6«E#(jn + 3i:)(i-io) • 



(5.4b) 



(5.5a) 



(5.5b) 



n=0 



B. Vector perturbations 



For vector perturbations ( |2.14| ) reads: 



x- 



XX 



^(x') + E7i 1) (jo+J 2 ) (n) (x-io 



n=0 



+ (1 - «H P f . (5.6) 



The covariant conservation of the energy-momentum tensor ([2.15|) is valid for the gas and 
the fluid seperately, i.e., v c p f = const. In terms of (5^ the infinite sum reads: 



3 00 

^H^n Un+j'n) (x - X ) 

Z 71=0 
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C. Tensor perturbations 



Since the anisotropic pressure for the perfect fluid component vanishes, ( |2.16| ) yields 



x 2 H" + -H' + H 



x 



-24a 



X dx'f X - X jH'(x- 

xo (X — X 



A2 



00 1 2 1 



(5.7) 



In terms of flffl the inhomogeneity reads: 

00 

3«E^ 2) (^ + 2 ^+^ y) 



n=0 



yX - X Q ) 



(5.8) 



The relation between the coefficients (3^ and 

7 (a) is 

the same as in the purely collisionless 

case. 



VI. SOLUTIONS 

The above integro-differential equations cannot, in general, be solved by, say, a power 
series ansatz about xo- Instead various methods have been employed to solve them approx- 
imately: direct numerical integration [|TT|.|T^,^]; approximation by a finite sum of spherical 
Bessel functions |L0) ; or by a finite system of ordinary differential equations . 

However, for x — > 0, which usually is the most interesting point to define initial con- 
ditions anyway, a generalized power series ansatz can be solved recursively. Exact regular 
solutions were found in this way in Refs. ||T6 . 19|] , where it turned out that the power series 
involved have infinite radius of convergence, and converge faster than trigonometric func- 
tions do. In the following, this will be generalized to include singular solutions as well as 
the most general initial data at Xq = 0. Some of the more unwieldy details are relegated to 
the Appendix. 

If F(x) represents one of the metric potentials $, II, H, the general solution with 
x = turns out to be of the form 



F(x) = dF t 



rcg 



X 



C 2 x u F t 



sing 



X 



with 



reg,smg 



n=0 



(x) = £ c r «^- 
nl 



(6.1a) 



(6.1b) 



and Ci, C2 arbitrary constants. This form follows from demanding that the coefficients (5^ 
remain finite as Xq — ► 0. 

When inserted into the integro-differential equations ( |5.4| . |5.6| . [5~7| ) this leads to solvable 
recursion relations for the coefficients c n such that c n = c n (co, . . . , c n -\\ a, {/?}), because the 
power series representation of the spherical Bessel functions, 
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i\2m+k 



7 , ( x — x ') = : I : L (Q 2) 

^ (2m)!(2m + l)(2m + 3)---(2m + 2£; + l)' v ' ; 



leads to integrals 



HI I 

dx'fa; - x')V = 7 ; -r 1 rx 1+m+u , (6.3) 

where m is a natural number and i/ arbitrary. 

In solving these recursion relations it turns out that the coefficients (3^ (or, equivalently, 
7^) in Eqs. ( |5.4y5.6| , |0| ) either have to be all zero, which also puts to zero F reg except for 
H ieg , or have to satisfy certain constraints dictated by the values F reg (0). In each case a sin- 
gular part in Eq. ( |6.1a| ) exists that can be superimposed on the regular solutions by choosing 
C2 ^ 0. With xo = 0, only the regular part depends on the initial data. The presence of 
singular perturbations means that the initial singularity is no longer approximately one of 
FLRW type, but is essentially anisotropic. Linear perturbation theory then applies only for 
those values of x where F has become sufficiently small. 

It turns out that the singular behavior for x — > is given by 

(6.4) 

with a = 0, 1, 2 for (a) = S, V, T. If the fraction a of energy density contained in collisionless 
matter exceeds a crit . = ^, o becomes complex, leading to an essential singularity at x = 
with F oscillating like cos([Ima] lnx). This asymptotic behavior has been found previously 
by Zakharov [12] and examined further by Vishniac []13|.Q As explained by the latter, the 




essentially singular, oscillatory behavior for x — > arises because the momentum flux carried 
by the collisionless particles depends preferentially on the expansion of the universe in their 
direction of travel, instead of the net expansion of the volume, which leads to an excessive 
feedback for a > a C rit.- Thus in this case, the initial singularity is of a mixmaster type 
33 1 . It has been shown by Lukash et al. []3Jj| that collisionless matter provides an efficient 



means to disperse the initially strong anisotropy of such models; the singular parts of our 
solutions (|6.1a|) could therefore correspond to the later stage of such a scenario where the 
linear regime has finally been reached. 



A. Scalar perturbations 

Regular solutions for scalar perturbations, where <&(x) and U(x) stay finite for x — > 0, 
are possible only when the coefficients (3^ in Eqs. ( [5.4| ) and (|5.5|) satisfy certain constraints. 
Evaluation of Eqs. (|5.4|) and their first derivatives at X — Xq = yields restrictions on the 
initial matter distribution, i.e., 



1 The result of Eq. (|6.4|) agrees with the asymptotic behavior found in Ref. [12| for scalar, vector, 



and tensor perturbations, but not with Ref. [13], where the vector case was reported to have a 
different critical value a cr it.. 
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Pi 0) = 4 0) = 4 0) = o , 



(6.5) 



and relations between the initial values of the metric potentials 

$(0) + 211(0) = -4/^ 0) (6.6a) 

and 

2$'(0) = -n'(O) . (6.6b) 

The restrictions on (3^\ /?{ , and (3^ can be understand immediately from the Einstein 
equations written in the form of Eqs. ( [2.11b| ) and ( [2.12b| ) together with covariant conserva- 
tion of the energy- momentum tensor, Eq. ( [2.13| ), by inserting the definitions in Eqs. Q3.9| ) 
and ( p. 13 ); the vanishing of (3^ is a consequence of demanding regularity of II in Eq. ( 5.4b|) . 



The explicit recursion relations defining $ c n, $ p f , and II follow from the Eqs. ( ^.12a| ) and 



( [5.4| ). They are listed in the first part of the Appendix. The recursion relations determine 
the regular parts of all potentials uniquely for a given set of respecting Eq. ( |6.5|) (or 

equivalently a set of 7^°' ) 's, satisfying more involved conditions). The singular solutions do 
not depend on the inhomogeneous terms. Thus they can be added to any regular solution. 
Singular solutions for the metric potentials with poles of first or second order are forbidden 
for xq = since the integral (|6.3|) is then divergent for all negative integers v. For finite xq, 
however such solutions exist. (An example can be found in Ref. Eq. (40b).) 



In Fig. 1, regular solutions for e m , which coincides with the energy- density contrast S on 
comoving hypersurfaces, are given in the purely collisionless case for three different initial 
data: one where only f3^ is nonvanishing, a second where the first four coefficients have 
been put to zero, and a more generic one. All of them grow according to a power-law on 
superhorizon scales, where the second (dashed-line) solution exemplifies that an arbitrarily 
high power can be achieved by appropriately contrived initial date; on sub-horizon scales 
all solutions undergo damped oscillations whose amplitude falls off like 1/x eventually. The 
exact matching of the allowed asymptotic behaviors is however seen to depend strongly on 
the form of the initial data. The asymptotic regimes themselves are determined by rather 
different physical situations. On superhorizon scales, i < 1, everything can be viewed 
as being determined by the global geometry set up in accordance with the inital matter 
distribution. For i> 1, microphysics becomes important, and the damped oscillations there 
can be understood as the dispersion of density contrast carried by the collisionless particles 
through directional dispersion |14| . The intermediate growth of the amplitude of the dotted 



solution in Fig. 1 after its 5th maximum is due to a nonzero (3$ . This demonstrates that 
higher momenta in the initial distribution function may have considerable effects even at 
subhorizon scales. 

In Fig. 2, a particular regular solution is shown in the two-component case with equal 
energy density in the perfect-fluid and in the collisionless component. In order that the 
collisionless component can be matched smoothly to a perfect-fluid behavior for x — > (in 



The sometimes used 14-moment approximation [25] therefore potentially misses important de- 
tails in the evolution of cosmological perturbations. 
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reality there will be a finite time and therefore a finite value of x 7^ 0, where the collisionless 
component decouples), the initial data have been chosen so that n(0) = 0, which makes 
^pf(O) = $cii(0). After horizon-crossing, the energy-density contrast that remains in the 
perfect-fluid component is seen to be undamped and to have a smaller phase velocity, namely 
c/y/3, in contrast to c in the strongly damped collisionless component. The ratio of 5 C \\ to 
5 p f at horizon crossing depends strongly on the initial data chosen. In Fig. 2 the collisionless 
component dominates over the perfect-fluid component up to fairly large values of x\ the 
solutions presented in Refs. |20| show the opposite behavior. [] This is due to having chosen 
/?5°' ) 7^ in the case of Fig. 2. 

In Figs. 3-5, singular solutions are displayed. These do not depend on the initial matter 
distribution at xq = 0, which underlines their geometrical nature. In Fig. 3, the energy- 
density contrast 6 and the anisotropic pressure ttt is plotted for the purely collisionless case 
a — 1. There are now oscillations also on superhorizon scales as discussed above, which 
go over to damped oscillations with phase velocity equal to c on subhorizon scales, showing 
the same asymptotic behavior there as did the regular solutions. In Fig. 4, the singular 
solutions are plotted for a = 5/32, where the perfect-fluid component is just strong enough 
to eliminate the superhorizon oscillations. For smaller values of a, there a two essentially 
different solutions with different degree of singularity. For a = 1/10, they are rendered in 
Fig. 5. In the two-component cases of Figs. 4 and 5, the perfect-fluid component turns out 
to have a regular limit for x — > 0, growing until horizon crossing, after which they again 
turn into undamped acoustic waves with phase velocity c/ V3. 



B. Vector perturbations 

The requirement of regularity yields the conditions 

^7 (/# > + l^P) = v c Pf (6.7a) 

and 

= (6.7b) 

for the vector perturbations. These conditions follow from Eq. ( |5.6j ) and its first derivative 
at x = xq = 0. The singular solutions may be added as in the scalar case. 

For the perfect-fluid case a = there are no regular solutions as mentioned at the end 
of Sec. II. However, with collisionless matter, vorticity can be generated on super-horizon 
scales, which after horizon-crossing dies out through directional dispersion. This is shown 
for various initial data in Fig. 6 for a = 1/2. The growth in v ~ x 2 at superhorizon scales 
is brought about by approaching a constant for x — » 0. A more rapid growth arises 
when \1/ itself goes like a positive power of x, which is the case for ^q 1 ^ = and nonvanishing 



3 The author of Ref. [26] claims that tt^ < e m in general, which is seen to hold true for certain 
initial conditions only. 
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higher coefficients. One such example is shown by the dashed line in Fig. 6. In all these 
cases there is a genuine production of vorticity, which in the perfect-fluid case is forbidden 
by the Helmholtz-Kelvin circulation theorem |!6] . 



With different initial matter distributions it is also possible to arrange for a nonvanishing 
initial vorticity in the perfect-fluid component, which is compensated on superhorizon scales 
by an equal amount with opposite sign in the collisionless component. A net vorticity then 
survives on subhorizon scales after the dispersion of the collisionless part. In Fig. 7 such 
solutions are plotted for several initial conditions. This interesting possibility has been 
studied extensively in Refs. JH||2Ti|], to which we refer the reader for more details. Here 



we just mention that the vorticity in the perfect-fluid component which by the presence of 
collisionless matter can be reconciled with an initial FLRW-type singularity, gives rise to 
magnetic fields during the transition to the matter-dominated area, and the limits set by 
the anisotropy of the cosmic microwave background are such that astrophysically interesting 
amounts of primordial magnetic fields seem possible . 

In Fig. 8, we display the singular vector solutions for a = 1, 5/32, and 1/10. For 
a > 5/32, there are again superhorizon oscillations, which cease at a = 5/32, and for 
smaller a two different power-law asymptotics occur. While v — > for x — > 0, the metric 
potential \I> in fact diverges. Hence, these solutions do not correspond to a Friedmann-type 
initial singularity, but are essentially anisotropic. 

On sub-horizon scales, where the mechanism of directional dispersion becomes operative, 
both the regular and the singular solutions decay like 1/x 2 for large x. 



C. Tensor perturbations 

A necessary condition for the regular solutions is obtained from Eq. (|5.7| ) at x = xq = 0, 

i.e., 

4 2) + ^ + ^? = • (6.8) 

In contrast to vector and scalar perturbations the initial value H(0) is not fixed by the 
initial matter distribution. This corresponds to the existence of gravitational waves. Both 
polarisations evolve with the same H(x), since there is no preferred direction. Again the 
singular solutions can be added to the regular ones. 

In Fig. 9 various regular solutions for the metric potential H are given. The ones that 
start out as constants at superhorizon scales are similar to the ones found already in the 
perfect-fluid case (see above), the only difference being that the amplitude drops more 
strongly at horizon-crossing. On subhorizon scales, these solutions become gravitational 
waves and their decay ~ 1/x is determined only be the expansion, exactly as in the perfect- 
fluid case. A novel type of solution is obtained, however, by putting (3q = 0. As shown by 
the dashed curve, this gives a growing solution on superhorizon scales. 

Finally, in Fig. 10 we exhibit also the singular solutions for H for the same set of param- 
eters as in Fig. 8. Again, the spectacular behavior is restricted to superhorizon scales, upon 
horizon crossing these solutions again describe ordinary, albeit primordial, gravitational 
waves. 
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VII. CONCLUSION 



We have presented a detailed analytical study of cosmological perturbations of an 
Einstein-de Sitter universe filled with collisionless massless particles and a perfect radia- 
tion fluid in terms of the gauge-invariant variables of Bardeen. A closed set of equations 
for the latter has been derived within the recently proposed approach based on thermal 
field theory on a curved background. There the response of matter under metric pertur- 
bations is determined by the thermal graviton two-point function, its leading terms in a 
high-temperature expansion describing ultrarelativistic collisionless matter. We have found 
complete equivalence with classical kinetic theory, when the latter is recast in terms of gauge 
invariant quantities along the lines of Ref . . A conceptual advantage of the thermal- field- 
theory approach appears to be its purely geometrical character — no explicit recourse to 
the generally gauge-variant distribution functions of kinetic theory is required. This might 
turn out to be useful when going beyond the case of collisionless matter, an issue that we 
intend to explore in a future work. 



In extension of the work of Ref. 16,19 , we have given the general solutions of the 



linear perturbation equations in terms of generalized power series expansions, the asymptotic 



behavior of which has been obtained previously in Ref. (T^]. In a strictly FLRW setting, 
only the regular solutions are admitted, but the general solutions might be of interest when 
to be matched to earlier, inflationary epochs, where the initial perturbations may have been 
generated quantum-mechanically []3"7[] . 

With special emphasis on the role of the initial data, we have also presented a selection of 
explicit results for regular and singular solutions for scalar, vector, and tensor cosmological 
perturbations. Typically, the regular solutions exhibit growth on superhorizon scales ~ x n , 
with n a characteristic integer, which can be increased by special choices of the initial data. 
The singular solutions have a small- x behavior dictated only by the ratio a which gives the 
ratio of the energy-density in collisionless matter over the total one. For a > 5/32, there is an 
essential singularity for x — > 0, giving rise to superhorizon oscillations, whereas for a < 5/32, 
all perturbations have a (singular) power-law behavior for small x. After horizon-crossing, 
all types of perturbations that are carried by collisionless matter decay. In the scalar and 
vector case, this is due to directional dispersion, whereas tensor perturbations only decrease 
according to the expansion of the FRLW universe. We have found that the transition between 
superhorizon and subhorizon regimes depends rather strongly on the initial data. Whereas 
previous investigations have always considered only the simplest cases, our results indicate 
that the effective matching of the asymptotic regimes can vary appreciably by allowing for 
more complicated initial matter distributions. 
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APPENDIX: 



To solve Eqs. ( |3.12a[ ), (ggj), ( gg ) and (^) the metric potentials $,$ c ii,n,*, and H 
are expanded into (generalized) power series. As mentioned in Sec. VI regular and singular 
solutions can be derived in such a way. Since the latter equations conserve the parity of 
the metric potentials, the regular solutions are calculated for the even and the odd part 
seperately. 

The recursion relations have been implemented in a Mathematica [38] code to calculate 
the first 60 coefficients and to plot the solutions in Figs. 1-10. 



1. Scalar perturbations 

The metric potentials are expanded into power series: 



$ = 


oo 

E( 

n=0 


-1 


™2n 

(2n) 


n = 


oo 

E( 

n=0 


-1 


™2n 

(2n) 


$cll = 


oo 

£( 

n=0 


-1 


™2n 

(2n) 



,. even 



+ 77r^rC dd + ^ (0) C ns , (Ala) 



(2n ) !V <Ven+ (^ <dd + ^° )< i ' (Mb) 

C en + (2^1) + ^%r s ) • (Ale) 

For simplicity e n := 2(0 n + 7r n ) (for regular even and odd solutions and for the singular 
solution) is defined. From Eq. ( |3.12a| ), the relations 

n-l 

3(2n + l)C cn = 30o vcn - 4 wn ~ 3er n - E ( 2/ + 4 ) e r n , (A2a) 

i=i 

n-l 

3(2n + 2)0° dd = 2(30° dd - e° dd ) - 3e° dd - ^ (21 + 5)e^ dd , (A2b) 



i=i 



3(cr(°) + 2n + l)0f S = ry+^tl) ^ " " 3e " nS " 

_ (2n)\ ^ r(cTW+2/ + l)(cTW + 2/ + 4) sing 

r(cr(o) + 2n + l) ^ (2/)! < 1 } 

determine with n > 1. For the odd regular solutions 2</>Q dd = — eg dd and the singular 
solutions obey (a^ + l)0o mg = — e o mg - Together with (n > 0) 

2/ + 2 

(2n + 2)(2n + lR ven = -6a £ T e£? + 



^ (2/ + 3)(2/ + 5) 
+6a(-l) n 



(2n + 2)!£ J,,^ fr 2 ^,, 4 



3(2n + 4)!^ 1 L )^n-2i+4 



, (2/)!!(4n-2/ + 5)!! 

/#)_.. 1 

(A3a) 



^ (2/)!!(4n-2/ + 9)!! 
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(2n + 3)(2n + 2)7r, 



odd 



n 91 4-9 



(2/ + 3)(2/ + 5) 



+6a(-l)' 



n+l 



-lV/3 (0) 



(2n + 3) 'S (2Z)!!(4n-2Z + 7)!! 



71+2 



+3(2n + 5)!^ 
(2n + 2)(2n + l)< ing = -6a 



-lV/3 {0) 



^ (2/)!!(4n-2/ + ll)!! 
(2n + 2)! " r(<r<°> +2n-2/ + l)(2/ + 2) 



E 



r(cr(°) + 2n + 3) ^ (2n - 2/)!(2/ + 3) (2/ + 5) 



(A3b) 
C! (A3c) 



from Eq. ( 5.4b|) , $ and II can be calculated. The condition (|6.5| ) has to be satisfied to gain 
regular solutions. The collisionless part of $, $ c n, is determined by Eq. ( |5.4aj) . From the 
constant term of the latter equation 

C en + 2 < cn = -4/ ? o° ) 
follows for the initial values. The higher terms give the relations (n > 0): 

1 3 \ „„_ 



(2n + 2)(2n + l)^ 



Qa 



i even even 

} n+l e n+l 



+ S V 2Z + 1 ' (2/ + 3)(2n + 3) 



■n-i+i 



+ 



+12a(-l) 7 



(2n + 2)!E 



2n-2i+2 



^ (2/)!!(4n-2/ + 5)!! 



-3(2n + 2)!(2n + 4)E 



n+2 (_1)^(°) 



2n-2Z+4 



(2n + 3)(2n + 2)<£ 



odd 



6a 



iodd pOdd 

■'n+l ~ e n+l 



^ (2l)\\(4n-2l + 9)\\ 
n+i / x 3 



+ E 



+ 



+12a(-l) r 



l= ^\2l + l (2/ + 3)(2n + 4) 

n+l 

[ h (2Z)!!(4n-2Z + 7)!l 

(-1)^ 



D odd 



(A4a) 



+ 



-3(2n + 3)!(2n + 5)E ' " "'' 



^ (2/)!!(4n-2Z + ll)!! 



(A4b) 



(2n + 2)(2n + l)gr s = 6a 



, sing sing 

e ra+l i 



(2n + 2)! 



■X 



r(a(°) + 2n + 3) 

3 \ r(o-(°) + 2n - 21 + 3) 



(A4c) 



X ^ (2/ + I + (2Z + 3)(a(°) + 2n + 3)y! (2n-2Z + 2)! e "- m _ 



sing 



The perfect fluid component is provided by $ p f = $ — $ 



'ell- 



2. Vector perturbations 



Accordingly, the metric potential \I/ is split up: 
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X 



2?i 



= V(-i) n - 



C ven + 



(2ra + 1) 



< dd + x CT(1 Vf g 



(A5) 



The regular even solutions are defined by 



(2n + 2)(2n + 1) + 



SO 



-24a £ 



even 
n—l 



(2Z + 3)(2Z + 5) 



+ 



+ 12a(-l) n 



n+l 



(2n + 2)!£ 



■l) l B (1) 

L ) H2U-21+2 



n+2 



^ (2Z)!!(4n-2Z + 5)!! 



-1V/3 (1) 

^ M2n-2Z+4 



+ (2™ + 4 ) ! E ( 2 /)!!(4n-2/ + 9)!! 
for n > 0. The solution for the regular odd part is provided for n > by 



(A6a) 



(2n + 3)(2n + 2) + ^ 
5 . 



^ dd = -24a]T 



odd 



J (2/ + 3)(2/ + 5) 



+ 



+ 12a(-l)' 



n+i (-lypP 



(2n + 3)!E 7 ^? 



2n-2«+3 



^ (2Z)!!(4n-2Z + 7)!! 



,,-■2 t_iypW 



+ (2- + 5)!£-A 



2n-2/+5 



^ (2Z)!!(4n-2Z + ll)!! 



(A6b) 



Additionaly the conditions (|S.7|) have to be satisfied. Singular solutions are fixed by 
and 



n{2n + 3) + 2na (1) J VC g = 

12 (2n)! A !>(!> + 2n - 2/ + 1) tff 

r(a( 1 ) + 2n + l) ^ (2n-2Z)! (2Z + 3)(2Z + 5) 



(A6c) 



for n > 1. 



3. Tensor perturbations 

For the metric potential if the ansatz 



x 



2n 



H = Y(-iy / , 



-even , 

v " (2n + l)'" n 



is made. Regular even solutions are defined by 

la 



(2n+ l)2n + 



n-l 



i even 
n n-l 



(A7) 



/£ ven = 2n(2n - l)C e i - 24a V 

V ' n ~ x fr[ (2Z + l)(2/ + 3)(2/ + 5) 



+ 



+ 3a(-l) n 



2n-2£ 



(2n) 'S (2Z)!!(4ra-2Z + l)!! 
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n+1 



L ) l J 2n-2l+2 



+ 2(2 " + 2)! S (2/)!!(4n-2/ + 5)!! 

71+2 



+ 



+ (2n + 4)!]T 



-l) l B {2) 

L ) l J 2n-2l+A 



^ (2Z)H(4n-2Z + 9)!! 



(A8a) 



for ri > 1. The initial value h^ en specifies the regular solution together with the inhomo- 
geneous terms. The latter have to respect the condition ( |6.8| ). The solution for the regular 
odd part is provided for n > by 



8a 



(2n + 2)(2n + 1) + — K Ad = (2n + l)2nh° n A \ - 24a £ 



Lodd 
""n-l 



+ 3a(-l) n 



(2n + l)!£ 



fr{ (2/ + l)(2/ + 3)(2/ + 5) 



+ 



J 2n-2l+l 



% (2Z)!!(4n-2Z + 3)!! 



n+l 



(-l)'/3 : (2) 



2n-2l+3 



2(2 " + 3) 'S (20!!(4„-2/ + 7)!! 



+ 



n+2 



+ (2n + 5)\J2 



-lV/3 (2) 



£ (2Z)!!(4n-2Z + ll)!! 



(A8b) 



Singular solutions are fixed by /ig mg and 



n(2n + 1) + 2na (2) ] /if s = n(2n - l)h^- 



-12a 



(2n)! 



E 



r>( 2 ) + 2n - 2Z + 1) 



LSing 



r(cx( 2 ) + 2n + 1) ^ (2n-20! (2/ + l)(2/ + 3)(2Z + 5) 



(A8c) 



for n > 1. 
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FIGURES 



FIG. 1. Regular density perturbations on comoving hypersurfaces in the purely collisionless 
case (a = 0). The full line shows a regular solution for \5\ = \e m \ with (3^ = —5/28 and all other 
Pn s vanishing. The dashed line represents a solution with [3^ = —1309/20 nonvanishing only. 
For the dotted line P^ takes the above value and P^ = —10, P%q = 60. 



FIG. 2. Regular scalar perturbations in a two-component universe (a = 1/2). The various lines 
show regular solutions for | <5 c il I (full line), |5 p f| (dotted line), and |vr^| (dashed line). The initial 
conditions are specified by /?J 0) = -l/4,/?( 0) = -1309/20, and pf ] = 7/8. 



FIG. 3. Singular solutions for |<5| (full line) and |7Tj. | (dashed line) for a = 1 (purely collisionless 
case). Superhorizon oscillations show up for x/tt < 1. Notice the logarithmic scale in i — the 
solutions are essentially singular as x — ► 0. The second solution corresponding to in Eq. fl6.4|) , 
which is not plotted, differs essentially by a phase. 



FIG. 4. Singular scalar perturbations for the critical value of a = 5/32. |<5 c u| (full line), \5 p f 
(dotted line), and |7Tj^| (dashed line) are plotted. 



FIG. 5. Singular scalar perturbations in the subcritical case (a < 5/32). The two different 
modes (<x±^) are shown by the full lines for \5 C \\\ and by the dotted ones for |<5 p f|. They have 
different power-law behavior on superhorizon scales. 



FIG. 6. Regular vector perturbations with v c p f = 0. The solutions for \v c \ are given 
with only = -21/16 nonzero (full line), p{ 1] = -1197/80 alone (dashed line), and 

pP = -21/16, p{ 1] = -10 and p{^ = 20 for the dotted line. 



FIG. 7. Regular vector perturbations with non-vanishing vorticity in the perfect-fluid compo- 
nent. The contribution of v c p f = —21/16 to \v c \ is plotted dashed-dotted. The other lines show 
various solutions for \v c \ using the same initial conditions as in Fig. || exept for P^\ which is 
determined by Eq. ( 6.7a| ). 



FIG. 8. Singular vector perturbations in terms of \v c \ . The different lines correspond to different 
values of a and to different modes. The solutions with, e.g., a = 1 (thin full line) show superhorizon 
oszillations, whereas solutions for a C rit. (thick full line) and below, e.g., a = 1/10 (dashed-dotted 
and dotted) do not. The latter two lines correspond to the two different modes from Eq. ( |6.4j ). In 
the supercritical case both modes differ by their phase only. 
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FIG. 9. Regular tensor perturbations. The solutions, representing gravitational waves for 

(2) 

x 1, are normalized to H(0) = 1, exept for the dashed line (JT(0) = 0,f3{ > = 63/2). The 

(2) 

full line shows the solution with all /% 's vanishing. Again the dotted line represents a more 

(2) (2) (2) (2) 

general case with /3\ = 1, ft\ = 5, /3g = —10, and /3y = 20 nonzero. 

FIG. 10. Singular tensor perturbations. \H\ is plotted for the same values of a as for the vector 
perturbations of Fig. 
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